Assignment of chromosomal locations for unassigned SNPs/scaffolds based on pair-wise linkage disequilibrium estimates

<p>Abstract</p> <p>Background</p> <p>Recent developments of high-density SNP chips across a number of species require accurate genetic maps. Despite rapid advances in genome sequence assembly and availability of a number of tools for creating genetic maps, the exact genome location for a number of SNPs from these SNP chips still remains unknown. We have developed a locus ordering procedure based on linkage disequilibrium (LODE) which provides estimation of the chromosomal positions of unaligned SNPs and scaffolds. It also provides an alternative means for verification of genetic maps. We exemplified LODE in cattle.</p> <p>Results</p> <p>The utility of the LODE procedure was demonstrated using data from 1,943 bulls genotyped for 73,569 SNPs across three different SNP chips. First, the utility of the procedure was tested by analysing the masked positions of 1,500 randomly-chosen SNPs with known locations (50 from each chromosome), representing three classes of minor allele frequencies (MAF), namely >0.05, 0.01<MAF ≤ 0.05 and 0.001<MAF ≤ 0.01. The efficiency (percentage of masked SNPs that could be assigned a location) was 96.7%, 30.6% and 2.0%; with an accuracy (the percentage of SNPs assigned correctly) of 99.9%, 98.9% and 33.3% in the three classes of MAF, respectively. The average precision for placement of the SNPs was 914, 3,137 and 6,853 kb, respectively. Secondly, 4,688 of 5,314 SNPs unpositioned in the Btau4.0 assembly were positioned using the LODE procedure. Based on these results, the positions of 485 unordered scaffolds were determined. The procedure was also used to validate the genome positions of 53,068 SNPs placed on Btau4.0 bovine assembly, resulting in identification of problem areas in the assembly. Finally, the accuracy of the LODE procedure was independently validated by comparative mapping on the hg18 human assembly.</p> <p>Conclusion</p> <p>The LODE procedure described in this study is an efficient and accurate method for positioning SNPs (MAF>0.05), for validating and checking the quality of a genome assembly, and offers a means for positioning of unordered scaffolds containing SNPs. The LODE procedure will be helpful in refining genome sequence assemblies, especially those being created from next-generation sequencing where high-throughput SNP discovery and genotyping platforms are integrated components of genome analysis.</p

Developing a community-based genetic nomenclature for anole lizards

Background: Comparative studies of amniotes have been hindered by a dearth of reptilian molecular sequences. With the genomic assembly of the green anole, Anolis carolinensis available, non-avian reptilian genes can now be compared to mammalian, avian, and amphibian homologs. Furthermore, with more than 350 extant species in the genus Anolis, anoles are an unparalleled example of tetrapod genetic diversity and divergence. As an important ecological, genetic and now genomic reference, it is imperative to develop a standardized Anolis gene nomenclature alongside associated vocabularies and other useful metrics. Results: Here we report the formation of the Anolis Gene Nomenclature Committee (AGNC) and propose a standardized evolutionary characterization code that will help researchers to define gene orthology and paralogy with tetrapod homologs, provide a system for naming novel genes in Anolis and other reptiles, furnish abbreviations to facilitate comparative studies among the Anolis species and related iguanid squamates, and classify the geographical origins of Anolis subpopulations. Conclusions: This report has been generated in close consultation with members of the Anolis and genomic research communities, and using public database resources including NCBI and Ensembl. Updates will continue to be regularly posted to new research community websites such as lizardbase. We anticipate that this standardized gene nomenclature will facilitate the accessibility of reptilian sequences for comparative studies among tetrapods and will further serve as a template for other communities in their sequencing and annotation initiatives.Organismic and Evolutionary BiologyOther Research Uni

A new look at the LTR retrotransposon content of the chicken genome

BACKGROUND: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact. RESULTS: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor. CONCLUSIONS: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size

Adenosine induces growth-cone turning of sensory neurons

The formation of appropriate connections between neurons and their specific targets is an essential step during development and repair of the nervous system. Growth cones are located at the leading edges of the growing neurites and respond to environmental cues in order to be guided to their final targets. Directional information can be coded by concentration gradients of substrate-bound or diffusible-guidance molecules. Here we show that concentration gradients of adenosine stimulate growth cones of sensory neurons (dorsal root ganglia) from chicken embryos to turn towards the adenosine source. This response is mediated by adenosine receptors. The subsequent signal transduction process involves cAMP. It may be speculated that the in vivo function of this response is concerned with the formation or the repair and regeneration of the peripheral nervous system

Avian Pathogenic Escherichia coli (APEC) Infection Alters Bone Marrow Transcriptome in Chickens

Avian pathogenic Escherichia coli (APEC) is a major cause of disease impacting animal health. The bone marrow is the reservoir of immature immune cells; however, it has not been examined to date for gene expression related to developmental changes (cell differentiation, maturation, programming) after APEC infection. Here, we study gene expression in the bone marrow between infected and non-infected animals, and between infected animals with mild (resistant) versus severe (susceptible) pathology, at two times post-infection. We sequenced 24 bone marrow RNA libraries generated from the six different treatment groups with four replicates each, and obtained an average of 22 million single-end, 100-bp reads per library. Genes were detected as differentially expressed (DE) between APEC treatments (mild pathology, severe pathology, and mock-challenged) at a given time point, or DE between 1 and 5 days post-infection (dpi) within the same treatment group. Results demonstrate that many immune cells, genes and related pathways are key contributors to the different responses to APEC infection between susceptible and resistant birds and between susceptible and non-challenged birds, at both times post-infection. In susceptible birds, lymphocyte differentiation, proliferation, and maturation were greatly impaired, while the innate and adaptive immune responses, including dendritic cells, monocytes and killer cell activity, TLR- and NOD-like receptor signaling, as well as T helper cells and many cytokine activities, were markedly enhanced. The resistant birds’ immune system, however, was similar to that of non-challenged birds. The DE genes in the immune cells and identified signaling models are representative of activation and resolution of infection in susceptible birds at both post-infection days. These novel results characterizing transcriptomic response to APEC infection reveal that there is combinatorial activity of multiple genes controlling myeloid cells, and B and T cell lymphopoiesis, as well as immune responses occurring in the bone marrow in these early stages of response to infection

Genetic variation exists for telomeric array organization within and among the genomes of normal, immortalized, and transformed chicken systems

This study investigated telomeric array organization of diverse chicken genotypes utilizing in vivo and in vitro cells having phenotypes with different proliferation potencies. Our experimental objective was to characterize the extent and nature of array variation present to explore the hypothesis that mega-telomeres are a universal and fixed feature of chicken genotypes. Four different genotypes were studied including normal (UCD 001, USDA-ADOL Line 0), immortalized (DF-1), and transformed (DT40) cells. Both cytogenetic and molecular approaches were utilized to develop an integrated view of telomeric array organization. It was determined that significant variation exists within and among chicken genotypes for chromosome-specific telomeric array organization and total genomic-telomeric sequence content. Although there was variation for mega-telomere number and distribution, two mega-telomere loci were in common among chicken genetic lines (GGA 9 and GGA W). The DF-1 cell line was discovered to maintain a complex derivative karyotype involving chromosome fusions in the homozygous and heterozygous condition. Also, the DF-1 cell line was found to contain the greatest amount of telomeric sequence per genome (17%) as compared to UCD 001 (5%) and DT40 (1.2%). The chicken is an excellent model for studying unique and universal features of vertebrate telomere biology, and characterization of the telomere length variation among genotypes will be useful in the exploration of mechanisms controlling telomere length maintenance in different cell types having unique phenotypes

Genome-wide analysis reveals the extent of EAV-HP integration in domestic chicken

Background: EAV-HP is an ancient retrovirus pre-dating Gallus speciation, which continues to circulate in modern chicken populations, and led to the emergence of avian leukosis virus subgroup J causing significant economic losses to the poultry industry. We mapped EAV-HP integration sites in Ethiopian village chickens, a Silkie, Taiwan Country chicken, red junglefowl Gallusgallus and several inbred experimental lines using whole-genome sequence data. Results: An average of 75.22 ± 9.52 integration sites per bird were identified, which collectively group into 279 intervals of which 5% are common to 90% of the genomes analysed and are suggestive of pre-domestication integration events. More than a third of intervals are specific to individual genomes, supporting active circulation of EAV-HP in modern chickens. Interval density is correlated with chromosome length (P<2.31−6), and 27 % of intervals are located within 5 kb of a transcript. Functional annotation clustering of genes reveals enrichment for immune-related functions (P<0.05). Conclusions: Our results illustrate a non-random distribution of EAV-HP in the genome, emphasising the importance it may have played in the adaptation of the species, and provide a platform from which to extend investigations on the co-evolutionary significance of endogenous retroviral genera with their hosts

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described